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Long-term variability in the X-ray emission of RX J0720.4-3125 
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Abstract. We detect a gradual, long-term change in the shape of the X-ray spectrum of the isolated neutron star RX J0720.4- 
3125, such that the spectrum of the source can no longer be described as a blackbody spectrum. The change is accompanied by 
an energy-dependent change in the pulse profile. If the X-ray emission is influenced by the magnetic field of the pulsar, these 
changes in spectral shape may point to precession of the neutron star. 
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1. Introduction 

RX J0720.4-3125 faaberl et alJll997l) belongs to the group of 
radio-quiet neutron stars discovered by ROSAT. As these ob- 
jects are characterised by soft X-ray emission from the sur- 
face, it was hoped that high resolution spectroscopy of those 
sources with the Chandra's Low Energy Grating Spectrometer 
(LETGS) or XMM-Newton' s, Reflection Grating spectrome- 
ter (RGS) would reveal line features from a thin atmosphere 
that could provide information on the surface gravity of the 
neutron star and hence on the equation of state of the neu- 
tron star matter. However, deep observatio ns with XM M- 
Newton and Chandra of RX J0720.4-3125 JPaerels etalTg OOl: 



iKaolan et alJ2003l) and RX J 1856.5-3754 (Bur witz et all200ll) 
have in that respect been somewhat of a disappointment, as 
the spectra are almost perfectly well described by a blackbody 
spectrum. The only unusual feature of the X-ray spectra is that 
the blackbody emission implies a radius too small for a neu- 
tron star and, moreover, underpredicts the observe d optical and 
UV fl ux jKaplan et alJl2003l: iBurwitz et alJEoOll iMotch etaD 
120031) . This implies that another cooler component is present 
and that the X-ray emission only comes from part of the sur- 
face. 

RX J0720.4-3125 has a low absorbing column density 
(Nfj ~ 10 20 cm~ 2 ). is at an approximate distance of w 300pc 
(Kapl an et all 12003). and shows sinusodial pulsations w ith a 
period of 8.39 s faaberl et alJll997llCropper et al.ll200ll) . The 
upper limit to the period derivative P < nx 10~ 13 ss _1 , with n 
'a few' (K aplan et alJl2002l) . implies a magnetic field strength 
B < 3 v /nl0 13 G and a characteristic age r c = 0.5P/P > 
(1.3/ti) Myr. As RX J0720.4-3125 was considered to be a per- 
fect blackbody spectrum with a temperature of fcTeB ~ 86 eV 



Table 1. XMM-Newton observations of RX J0720.4-3125 



Orbit 


Observation 


Date 


Julian date 


Duration (ksec) 


0078 


0124100101 


13-05-2000 


2451677.9 


65 


0175 


0132520301 


21-11-2000 


2451870.3 


30 


0533 


0156960201 


06-11-2002 


2452585.4 


30 


0534 


0156960401 


08-11-2002 


2452498.5 


32 


0622 


0158360201 


02-05-2003 


2452762.5 


81 


0711 


0161960201 


27-10-2003 


2452940.5 


45 
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it was observed several times by XMM-Newton for calibration 
purposes. However, as we will demonstrate here, the spectrum 
of RX J0720.4-3125 has slowly hardened and cannot be de- 
scribed any more by a blackbody spectrum, while the pulse 
shape has become narrower and the phase dependence of the 

spectrum changed with time. 

While we were finalising this paper, lHaberl et alJ (|2003b) 
released a pre-print in which they show that the spectrum of 
RX J0720.4-3125 depends on the pulse phase. They also dis- 
cuss changes in the spectrum as found with EPIC, but assign 
these to calibration inaccuracies. 



2. Observations 

Table ^ gives the log of XMM-Newton observations of 
RX J0720.4-3125. Since the source spectrum is rather soft, and 
most of the emission comes from below ~ 1.8 keV, we use 
RGS ( den Her der et alJ2 001 ) data for the spectral analysis; due 
to the nature of the gratings, the effective area of the RGS is 
rather insensitive to changes in CCD gain and charge transfer 
inefficiency. The RGS was operated in normal spectroscopic 
mode, yielding a time resolution of 4 s, making it unsuitable 
for a prop er timing analysis. T herefore, for this purpose we use 
EPIC/PN (Striider et al. 2001) since it has the largest effective 
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area among the instruments on board XMM-Newton, and oper- 
ates in a high-time resolution mode. 

During the first four observations, PN was operated in full- 
frame mode, which produces an image of 378 x 384 pixels 
with a time resolution of 73.4 ms. During the last two obser- 
vations the PN was switched to small-window mode, yielding 
images of 63 x 64 pixels with a time resolution of 6 ms. The 
thin filter was used in revolutions 0078, 0533, and 0534, the 
medium one in revolution 0175, and the thick one in revolution 
0622. The observation of revolution 071 1 started with the thin 
filter, and later switched to the medium one. All data were pro- 
cessed using the standard XMM-Newton software S AS version 
5.4.1. 

2.1. Spectral Analysis 

Figure[0shows the RGS spectra for the 6 revolutions. The in- 
divual RGS's agree with each other within the statistical un- 
certainties, so the data from both RGS's have been combined. 
From the Figure it is apparent that the spectrum of RX J0720.4- 
3125 gradually changes, while the total flux appears to remain 
fairly constant. The change of the spectrum is most noticeable 
in orbit 07 1 1 , where there is an increase of the flux in the 1 5 A 
to 25 A range, and a drop in the flux above 30 A. 

To verify that the RGS instrument is stable and that its ef- 
fective area does not change, we have also analyzed data from 
two sources that are stable on a time scale of years. 

The first source is the supernova remnant IE 0102-72.3 in 
the Small Magellanic Cloud, a calibration source that is regu- 
larly observed with the RGS. From the strength of the bright 
emission lines from this source in the range 13 to 23 A, as 
measured in data from orbits 0065 to 071 1, we find that the ef- 
fective area of the instrument in this range changes by less than 
3%. In the long wavelength range the SNR 0102-72.3 spectrum 
has a weak carbon line at 33.7 A. The upper limit on possi- 
ble changes in the effective area of the RGS instrument at this 
wavelength is about 10%, mainly determined by the limited 
statistics in this line. All data from IE 0102-72.3 are consistent 
with no change in the RGS effective area. The second source is 
the calibration source Mrk 421. An upper limit to any change 
in the neutral Oxygen absorption edge around 22.8 A, which is 
a tracer of possible instrument contamination (e.g. ice), is 5%. 

Since all upper limits on possible changes in effective area 
are well below the changes seen in RX J0720.4-3125 we con- 
clude that the X-ray spectrum of RX J0720.4-3 125 itself is sub- 
ject to change over the course of about 3 years. 

For the s pectral analysi s we use the XSPEC spectral fit- 
ting package (Arnaud 1996). As noted bv lPaerels et alJ (12001 ) 
the data of orbit 0078 can be well fitted by an absorbed black- 
body spectrum. With the improved calibration of the RGS we 
find a temperature of fcTeB = 86.7 ± 0.3 eV and an inter- 
stellar column density of A H = (1.41 ± 0.07) x 10 20 cm" 2 , 
in reasonable agreement with the values obtained with the 
Chandra LETGS instrument, kTm = 81.4 ± 1.3 eV and 
A H = fl.32zfc0.14) x 10 20 cm" 2 , by Kap lan et aD ( l2003l) . The 
LETGS observation was made in February 2000, four months 
before the first XMM-Newton observation. 
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Fig. 1. Fluxed RGS spectra of the different observations of 
RX J0720.4-3125. The spectra are shown as the light grey 
histograms. The black lines are the smoothed data from orbit 
0078. The differences between the lines and the different data 
show the variations of the spectrum with respect to the data of 
orbit 0078. Especially the large change at orbit 07 1 1 with an in- 
crease of flux between 15 and 25 A can be clearly recognized. 

Fits with the same model to the remaining spectra yield 
a temperature increase with time, as well as changes as large 
as AAh ~ 5 x 10 20 cm~ 2 in the interstellar absorption. The 
model, however, does not provide satisfactory fits to all the 
spectra. 

Next, we fit all the data simultaneously. Since large vari- 
ations of the interstellar absorption over such a short period 
are unlikely, we assume that Ah remains constant; we there- 
fore couple this parameter for all observations, while we let the 
blackbody temperature and normalization vary freely between 
observations (Table|2ji. Although these fits serve to emphasize 
the gradual hardening of the spectra, half of the spectra, espe- 
cially that of re volution 07 1 1, are p oorl y fit by this model. 

Following iHaberl et alJ J2003al) and Ivan Kerkwiikl 12004). 
next we fit the data with a model that consists of a blackbody 
affected by additional absorption by a broad gaussian line, as 
might be expected in the case of cyclotron absorption. In these 
fits, we constrain the parameters of the blackbody model to be 
the same for all the observations, whereas the parameters of the 
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Pulse phase 

Fig. 2. Pulse profile of RX J0720.4-3125 in the 0.1-1.2 keV 
band (panels 1 and 3 from top to bottom) and hardness-ratio 
light curve (panels 2 and 4) for revolutions 0078 and 071 1 using 
EPIC/PN. The 0.1-1.2 keV light curves are normalised to 1 
at the maximum. Phase is defined when the full-band light 
curve reaches its maximum. The solid line is the best sinusoidal 
fit to the full-band light curve of revolution 0078. The same 
sinusoidal fit is overplotted on top of the full-band light curve 
of the other revolution for comparison. For both observations 
the thin filter was used. 



line are allowed to vary independently. Although this model fits 
the data well, in most cases we find that the width of the line 
is larger than its central energy, or the central energy of the 
line lies outside the spectral range covered by the data. In those 
cases, the tail of the broad gaussian just serves to attenuate the 
emission of the source at long wavelengths. Furthermore, we 
find no clear trend in the central energy or the width of the line. 

We find that the ratio of the spectra of the last and the first 
observation is close to a power law. We therefore fit the data to 
an empirical model that consists of a blackbody multiplied by a 
power law E T , all affected by interstellar absorption. While the 
index of the multiplicative power law is allowed to change be- 
tween observations, for these fits we constrain the parameters 
of the blackbody and the interstellar absorption to be the same 
in all observations. While it is difficult to assign a physical in- 
terpretation to this model, it provides an acceptable description 
of the data in the RGS range (10-38 A, see Table 2), it has 
fewer parameters than the gaussian absorption model and, in 
addition, the index of the power law increases steadily over the 
course of the observations. 



2.2. Timing analysis 

Starting from the raw data, we first produce a list of calibrated 
events. To reduce pile-up, in the next step we select only single 
events as well as events that are not affected by some of the 
imperfections (bad columns, hot pixels, etc.) of the CCDs. We 
extract events within a 39 arcsec circle centred on the source. 
We barycenter these events using the SAS routine BARYCEN 
version 1.13.4, and we then separate the events according to 
their energy in 3 event lists; the bands that we use are 0. 1 to 1 .2 
keV, 0.1 to 0.4 keV, and 0.4 to 0.8 keV, respectively. 

For each observation we find the best period in the full band 
using an epoch folding technique; in all cases we find a pe- 
riod of 8.391 s , consistent w i th the value previously found for 
this source bv iKaplan et alJ J2003I) . We then produce folded 
light curves in the three bands, and we also compute a folded 
hardness-ratio light curve from the ratio of the 0.4-0.8 keV and 
the 0.1-0.4 keV light curves. In Figure|2]we show the 0.1-1.2 
keV and the hardness-ratio light curves. For each observation 
we define the phase such that the maximum of the full-band 
light curve occurs at phase zero; the phase of the hardness ratio 
light curves is the same as for the full-band light curves. 

The pulse profile in the 0.1-1.2 keV band, as well as the 
hardness-ratio pulse profile, change from one observation to the 
other. The first panel in Figure [2] shows a sinusoidal fit to the 
pulse profile during the first observation; the same sine function 
is overplotted to the full-band pulse profiles obtained from the 
other observations. It is apparent that the pulse profile becomes 
narrower with time. 

At the same time, the hardness-ratio pulse profile also 
changes. In the first observation there is a clear modulation, 
and the hardness-ratio profile leads the full-band light curve by 
0.064±0.017in phase. In the following observations the ampli- 
tude of the hardness-ratio modulation decreases and the phase 
difference between the full-band and the hardness-ratio light 
curves is consistent with zero. Eventually, in the last observa- 
tion the modulation increases again, but now the hardness-ratio 
light curve lags the full band-light curve by —0.126 ± 0.010 in 
phase. 

3. Conclusions and discussion 

The XMM-Newton data of RX J0720.4-3125 show that the 
spectrum of the source changes on a time scale of years, the 
first time ever that the X-ray spectrum of an isolated neutron 
star, other then soft gamma-ray repeaters or anomalous X-ray 
pulsars, is seen to change. Whereas the changes are most pro- 
nounced in the last observation, we think that the actual change 
is gradual, as witnessed by a gradual increase in the tempera- 
tures derived from the blackbody fits; or by a gradual increase 
in the index of the powerlaw in the fits with a blackbody mul- 
tiplied with a power law (Tableland Figure^. The spectral 
changes are accompanied by an energy-dependent change in 
the pulse shape; in particular the pulse phase where the spec- 
trum is hardest has moved with respect to the phase of maxi- 
mum flux (Figure|2j. 

The phase (i.e. angle) dependent spectrum of single neu- 
tron stars is currently not explained. The broad absorption 
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Table 2. Spectral fits a to the observations of RX J0720.4-3125 



Orbit 


kT (eV) 




pc 




Flux d 












(10-23 A) 


(23-38 A) 


(10-38 A) 


0078 


81.3(3) 


2.54(3) 


1.44(5) 


1.86(4) 


4.72(7) 


6.58(10) 


0175 


81.4(5) 


2.42(3) 


1.45(5) 


1.77(5) 


4.49(7) 


6.26(11) 


0533 


84.6(4) 


2.22(3) 


1.74(5) 


1.89(4) 


4.20(7) 


6.09(11) 


0534 


84.0(4) 


2.28(3) 


1.67(5) 


1.88(4) 


4.31(7) 


6.19(11) 


0622 


86.3(3) 


2.07(2) 


1.85(6) 


1.93(4) 


4.08(7) 


6.01(06) 


0711 


98.2(4) 


1.96(2) 


2.70(8) 


2.81(4) 


4.01(6) 


6.82(09) 



a Fit range 10-38 A. For the interstellar absorption component we use cross sections from Verner et al. > 1996), and abundances from Wilms et al. > 200C). Numbers in parentheses 
are 1-cr confidence limits in the last digit(s). For the fits with a blackbody, the best-fit value for -/Vh, constrained to be the same in all observations, is (4.3 ± 0. 1) X 10 20 cm — 2 . 
For the fits with a blackbody times a power law, the best-fit values for JVh , kT, and the bolometric luminosity of the blackbody, constrained to be the same in all observations are, 
respectively, (1.28 ±0.1) x 10 2Q cm~ 2 , 70.3 ±0.1 eV, and (3.04 ± 0.02) x 10 32 t2 2 00 erg s~\ with d 3 oo the distance to RX J0720.4-3 125 in units of 300 pc. For the 
fits with a blackbody model, \ 2 = 874.1 for 662 de grees of freedom. For the fits with a blackbody times a power-law model — 776.8 for 668 degrees of freedom. 

b Bolometric luminosity of the blackbody in units of 10 32 d^ 0Q erg s — 1 , where (L^qq is the distance to the source in units of 300 pc. 

c Index of the power law in the model that consists of a blackbody times a power law (see text) 

d Observed flux, obtained from the fits with a blackbody times a power-law; only statistical errors are indicated. The systematic errors in the flux are of order 5-10%. 



features have been interpreted as a proton-cyclotron absorp- 
tion feature dHaberl et all2003ah . In pulsars with a strong field 
(probably stronger than the limit for RX J0720.4-3125) the ab- 
sorption feature may be weakened by the strong-field quantum 
electrodynamics effect of vacuum resonance mode conversion 
jLai & Hol2003l) . The neutron star spectra have also been inter- 
preted as due to cyclotron-resonance scattering of the spectrum 
from the surface of the neutron star, by e lectron-positron pairs 
in the magnetosphere dRudermanll2003h . For both interpreta- 
tions the spectrum is likely to be angle and energy dependent, 
in accordance with the variation of the X-ray spectrum (as mea- 
sured by hardness ratio) with pulse phase. 

To explain the gradual, long-term variation we consider two 
general possibilities: either the intrinsic spectrum of the neu- 
tron star changes, or our view of the neutron star changes. The 
intrinsic spectrum of the neutron star could change, due to en- 
ergy release deep in the neutron star, due to a glitch for exam- 
ple, causing the surface to gradually become hotter. We con- 
sider this unlikely as the explanation for the changing spectrum 
of RX J0720.4-3125, because it does not explain the change in 
pulse shape. Also this model would predict the total flux to in- 
crease with the temperature, in contrast to what is observed. 
Another possibility, valid for Ruderman's model, would be 
that the electron-positron plasma surrounding the neutron star 
changes. So far, there is no specific prediction in the Ruderman 
model for changes in the magnetospheric plasma on a year-long 
time scale. 

Therefore, we suggest that the variation in the spectrum of 
RX J0720.4-3125 is caused by a change in the angle under 
which we see the emitting region and/or the covering electron- 
positron plasma, caused by precession of the neutron star. 
Precession arises when the form of the neutron star deviates 
from a perfect sphere and its rotation is not around a princi- 
pal axis (as reviewed by e.g. Link 2003). The changes in pulse 
form and phase of the radio pulsar B 1828— 1 1 are successfully 
described w ith fre e precession, and have a time scale of years 
(Stai rs et alJl20 00). i.e. comparable to the time scale on which 
the variations in the spectrum and pulses of RX J0720.4-3125 
occur. Whether a change in viewing angle can produce the ob- 
served hardening of the spectrum, with only a modest change in 



overall flux, will depend on the details of the emission model, 
and requires further investigation. 

To further investigate the cause of the long-term variabil- 
ity it will be useful to investigate the phase-dependence of the 
spectra in more detail; to obtain a reliable period derivative of 
the pulses, and look for the phase changes expected from pre- 
cession; and finally to see whether the changes continue. 
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